Electronic structure and band parameters for ZnX (X=0, S, Se, Te) 
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First-principles density-functional calculations have been performed for zinc monochalcogenides 
with zinc-blende- and wurtzite-type structures. It is shown that the local-density approximation 
underestimates the band gap, misplaces the energy levels of the Zn-3d states, and overestimates the 
crystal-field splitting energy. Without spin-orbit coupling, the order of the states at the top of VB 
is found to be normal for all the ZnX phases considered. Upon inclusion of the spin-orbit coupling 
in calculations, ZnO in zinc-blende- and wurtzite-type phases become anomalous. It is shown that 
the Zn-3d electrons are responsible for the anomalous order. The effective masses of electrons and 
holes have been calculated and found that holes are much anisotropic and heavier than the electrons 
in agreement with experimental findings. The typical errors in calculated band gaps and related 
parameters originate from strong Coulomb correlations, which are found to be highly significant 
in ZnO. The LDA+C/ approach is found to correct the strong correlation of the Zn-3d electrons, 
and thus improves the agreement with the experimentally established location of the Zn-3d levels. 
Consequently, it increases significantly the parameters underestimated in the pure LDA calculations. 
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I. INTRODUCTION 

Investigation of properties of zinc monochalcogenides 
[ZnX (X=0, S, Se, Te) with zinc-blende- (z-) and 
wurtzite-(w-)type structure] has promoted much interest 
because of their numerous applications in optoelectronic 
devices such as visual displays, high-density optical mem- 
ories, transparent conductors, solid-state laser devices, 
photodetectors, solar cells etc. Despite frequent studies 
over more than four decades many fundamental questions 
still remain open. One of them is the eigenvalue prob- 
lem related to a severe underestimation of the band gaps 
(Eg), energy levels of the Zn-3d electrons (Ed) (for ZnO 
0, ZnS U and ZnSe Q), splitting energies of the states 
at the top of the valence band (VB), spin-orbit (SO) cou- 
pling and crystal-field (CF) splitting energies (Aso and 
A CF , respectively) calculated according to the density- 
functional theory (DFT) within the local-density approx- 
imation (LDA). These problems have been the subject of 
numerous studies with different methods such as the LDA 
plus self- interaction correction (LDA+5/C) 0. [J. 
LDA plus the multiorbital mean-field Hubbard potential 
(LDA+C7) |(| (which includes the on-site Coulomb inter- 
action in the LDA Hamiltonian) , and different version of 
GW approximation 0, 0, 0- In the latter approxima- 
tion "G" stands for one-particle Green's function as de- 
rived from many-body perturbation theory and "W" for 
Coulomb screened interactions. This approach takes into 
account both non-locality and energy-dependent features 
of correlations in many-body system. However, none of 
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the above approaches has been able to remedy the eigen- 
value problem except the combination of exact-exchange 
DFT calculations in the optimizcd-cffective-potential ap- 
proach with GW , which is found to give better agree- 
ment with the experimental band gaps and the location 
of the Zn-3d energy levels. 

The order of the states at the top of the VB in ZnO-w 
is one of the topics which is still frequently debated 0,0- 
Another debated aspect is the effective masses of the 
charge carriers although these are more indefinite param- 
eters for the ZnX phases. At present the effective masses 
from different ab initio packages and experiments scat- 
ter appreciably in publications on ZnO 0, IU El El E3 
and ZnTe [l3|. In this work ZnX (X=0, S, Se, Te) in z- 
and w-type structural arrangements is studied by first- 
principles calculations within the LDA, generalized gra- 
dient approximation (GGA), and LDA+C/ approaches. 
See Ref. [HI for a more detailed account of our findings. 

II. COMPUTATIONAL DETAILS 

The electronic band structure of the ZnX phases is 
studied using the VASP-PAW package El> which cal- 
culates the Kohn-Sham eigenvalues by the DFT 09 
within the LDA GGA [3, and LDA+C/ [H |c| 
approximations. The exchange and correlation energy 
per electron have been described by the Perdew-Zunger 
parametrization [2l| of the quantum Monte Carlo pro- 
cedure of Ceperley-Alder [231 . The interaction between 
electrons and atomic cores is described by means of non- 
norm-conserving pseudopotentials implemented in the 
VASP package, which are generated in accordance to the 
projector-augmented- wave (PAW) method |2:il'2 1|. Self- 
consistent calculations were performed using a 10 x 10 x 
10 mesh frame according to Monkhorst-Pack scheme for 



z-type phases and to T-centered grids for w-type ] 
The completely filled semicore-Zn-3d states hav 
considered as valence states. Values for Acf, Ag 
the average band gap Eq are calculated within the 
cubic model |25l |26| . For band-structure calculati 
used the experimentally determined crystal-struct 
rameters for all ZnX phases considered. 

Without SO coupling the top of the VB for phas< 
w-type structure is split into a doublet T 5 and a sin 
state by the crystal-field. Inclusion of SO couplin 
rise to three twofold degenerate valence bands, wh 
denoted as A, B, and C states. The symmetry of 
these three bands are of V-? character and one of T 
acter. Band gaps E g , Ea, Eb, and Ec are defined 
difference between the conduction band (CB) mi: 
and energy levels of the A, B, and C states, respei 
Without SO coupling the VB spectrum near the 1 
for the ZnX-z phases originates from the sixfold 
erate state. The SO interaction splits the T] 
into a fourfold degenerate T s and doubly degenei 
levels. 

III. RESULTS AND DISCUSSION 

A. Value of the parameters U and J 

The simplified rotationally invariant LDA+U ap- 
proach |nj used in this work requires knowledge of the 
values of the parameters U and J. Since these param- 
eters do not explicitly take into account the final state, 
values of U and J were found empirically from LDA+i7 
band structure calculations as a function of U and J such 
that the value of Ed obtained for particular U and J fit 
with the experimentally determined location. For com- 
parison, the values of U and J have been calculated for 
some of the compounds within the constrain DFT [36| . 
showing that the calculated values to some extent agrees 
with those extracted semiempirically. 

B. Eigenvalues 

The results of the band structure calculations are 
listed in Table [I] It is found that the band gaps 
(Eg, Ea,Eb,Ec) and the mean energy level Ed severely 
underestimated in the LDA calculations. The error in the 
LDA calculated energies is quite pronounced for ZnO (see 
e.g. Fig.^for ZnO-z) compared to the other ZnA phases 
and the discrepancy exceeds the usual error for LDA cal- 
culations. Furthermore, compared to experimental data 
the CF splitting energy (Acf) is severely overestimated 
for ZnO-w by 2.4 times and underestimated for ZnS-w 
by around 1.2 times. Ago is overestimated for ZnO-w to 
about 12.2 times, underestimated for ZnS-w by about 3.2 
times, but agrees well with experimental data for ZnSe-z 
and ZnTc-z. Furthermore, Ago is found to be negative 
for ZnO-z and -w phases in agreement with results of 
Refs. HE and Ea IM El while it is positive for the 
other ZnA" considered. 

The GGA approach corrected the above mentioned 
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FIG. 1: Band dispersion for ZnO-z calculated within LDA 
(solid lines) and LDA+U (dotted lines) approaches. The 
Fermi level is set at zero of energy. 

LDA deficiency only to a little extent. However, LDA+U 
significantly increased the values of E g , Ea, Eb, and Ec, 
and shifted the energy levels of the Zn-3ci electrons up 
to experimentally determined limits, but only slightly 
changed the values of Ago and Acf for the ZnX phases 
except ZnO-z and -w. The dependence of Acf, Ago 
on U for ZnO-w is displayed in Fig. |2ta) and of Ago 
for ZnO-z in Fig. 13b). Analysis of these illustrations 
shows that Ago remains negative for U < 9.0 eV for 
ZnO-w and for U < 8.0 eV for ZnO-z. For higher values 
of U, Ago becomes a complex number for ZnO-w with 
a non-zero imaginary part (in itself physically meaning- 
less), and changes sign for ZnO-z from negative to posi- 
tive. 

The SO splitting energy is much smaller than 1.0 eV 
for all phases except ZnSe-z (Ago = 0.4 eV) and ZnTc- 
z (0.97 eV). Regardless of the computational approach 
used, the numerical values of Ago for all ZnX phases re- 
mained almost unchanged. The present values for Ago 
for the ZnX-z phases are in good agreement with theoret- 
ical calculations Ei| by the LAPW method and available 
experimental data. 

The variation of the energy splitting between the A, B, 
and C states, viz. Ea — Eb, Ea — Ec, and Ea — Ec, at 
the top of the VB of ZnO-w is studied as a function of 
U as shown in Fig. |U^c). For U < 9.0 eV the energy 
splitting decreases with increasing U. At higher values 
of U, Ea — Eb becomes negative and decreases, Eb — Ec 
increases, while Ea — Ec stays more or less constant. 

C. Density of states 

The density of states corresponding to the Zn-3d lev- 
els according to the LDA calculations are inappropriately 
close to the CB (see e.g. Fig.[3]for ZnO-z), which contra- 
dicts the findings from XPS, XES, and UPS experiments 
[13, IH, Hj| . Furthermore, these states and the top of the 
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TABLE I: Band gaps E g , E g (A), E g (B), E g (C), and Eq, crystal-field (A^p, Acf), and spin-orbit (Aso) splitting energies (all 
in eV) for ZnX phases with z- and w-type structure calculated within LDA, GGA, and LDA+C7 approaches. E g and A^p refer 
to calculations whithout the SO coupling, in all other calculations the SO interactions are accounted for. Experimental values 
are quoted when available. 



Phase 


Method 


-E 3 


Eg(A) 


Eg(B) 


E g (C) 


E 


E d 


A 

CF 


A C F 


Aso 


ZnO-w 


LDA 


0.7442 


0.7238 


0.7561 


0.8385 


0.773 


~ 5.00 


0.0951 


0.0928 


-0.0429 




GGA 


0.8044 


0.7831 


0.8165 


0.9000 


0.833 


~ 5.00 


0.0967 


0.0944 


-0.0443 




LDA+C/ 


1.9884 


2.0080 


2.0528 


2.0530 


2.038 


~ 10.00 










Expt. [26] 




3.4410 


3.4434 


3.4817 


3.455 






0.0394 


-0.0035 


ZnS-w 


LDA 


1.9896 


1.9681 


1.9947 


2.0734 


2.012 


~ 6.50 


0.0688 


0.0518 


0.0269 




GGA 


2.2322 


2.2114 


2.2361 


2.3104 


2.253 


~ 6.00 


0.0661 


0.0494 


0.0249 




LDA+C/ 


2.2828 


2.2595 


2.2858 


2.3662 


2.304 


~ 8.20 


0.0588 


0.0549 


0.0256 




Expt. [43] 




3.8643 


3.8932 


3.9808 








0.0580 


0.0860 




Expt. 




3.8715 


3.8998 










0.0055 


0.0920 


ZnSc-w 


LDA 


1.0704 


0.9389 


1.0080 


1.3789 


1.109 


~ 6.50 


0.1138 


0.3243 


0.0467 




GGA 


1.3271 


1.2004 


1.2678 


1.6243 


1.364 


~ 6.50 


0.1122 


0.3105 


0.0460 




LDA+C/ 


1.4039 


1.2708 


1.3336 


1.7209 


1.442 


~ 9.30 


0.1010 


0.3465 


0.0408 




Expt. ^,^3] 




2.8600 


2.8760 


2.9260 




~ 9.20 








ZnTe-w 


LDA 


1.0519 


0.7600 


0.8200 


1.6912 


1.091 


~ 7.50 


0.0857 


0.8383 


0.0332 




GGA 


1.2577 


0.9743 


1.0319 


1.8754 


1.294 


~ 7.20 


0.0835 


0.8115 


0.0320 




LDA+C/ 


1.2826 


0.9897 


1.0434 


1.8818 


1.305 


~ 9.50 


0.0754 


0.8088 


0.0296 


ZnO-z 


LDA 


0.5732 


0.5552 




0.5883 


0.577 


~ 4.60 






-0.0331 




GGA 


0.6409 


0.6152 




0.6494 


0.638 


~ 4.60 






-0.0342 




LDA+C/ 


1.4859 


1.4953 




1.4973 


1.496 


~ 7.90 






0.0020 


ZnS-z 


LDA 


1.8745 


1.8516 




1.9158 


1.873 


~ 6.10 






0.0642 




GGA 


2.1134 


2.0921 




2.1513 


2.112 


~ 6.00 






0.0592 




LDA+C/ 


2.3324 


2.3097 




2.3886 


2.336 


~ 9.00 






0.0789 




Expt. [431 




3.6800 




3.7400 




> 9.00 






0.0670 




Expt. ^ 




3.7800 




3.8500 












ZnSe-z 


LDA 


1.0793 


0.9484 




1.3409 


1.079 


~ 6.60 






0.3925 




GGA 


1.3349 


1.2089 




1.5858 


1.335 


~ 6.50 






0.3769 




LDA+C/ 


1.4214 


1.2908 




1.6995 


1.427 


~ 9.05 






0.4087 




Expt. [431 




2.7000 








~ 9.20 






0.4000 




Expt. ^ 




2.8200 














0.4003 


ZnTc-z 


LDA 


1.0607 


0.7715 




1.6681 


1.070 


~ 7.10 






0.8966 




GGA 


1.2671 


0.9862 




1.8533 


1.275 


~ 7.05 






0.8671 




LDA+C/ 


1.3287 


1.0456 




1.9561 


1.349 


~ 9.90 






0.9105 




Expt. [431 




2.3941 








~ 9.84 






0.9700 



Expt. ~ 10.30 



VB are hybridized. Distinct from the other ZnX phases 
considered, ZnO in both z- and w-type structure shows 
artificially widened Zn-3d states. These findings changed 
only slightly according to the GGA calculations. 

The LDA+C7 approach 0| was used to adjust the en- 
ergy levels of the Zn-3d electrons derived bands to exper- 
imentally established positions using semiempirical val- 
ues for the parameters U and J. Consequently, the band 
gaps calculated by this approach become more reason- 
able than the pure LDA-derived band gaps. Further, the 
height of the peaks in the DOS corresponding to the Zn- 
3d states calculated by the LDA+C/ becomes much larger 
than those calculated by the pure LDA. This indicates 
that the semicore Zn-3<i electrons have become more lo- 
calized than according to the pure LDA. For ZnO-z and 
-w the width of the Zn-3<i band calculated by LDA+C/ be- 
comes much narrower than that calculated by LDA. For 
the other ZnX phases, however, LDA+C/ only slightly 
changed the width of the Zn-3<i bands, which leads one 
to conclude that the Coulomb correlation effects for ZnO 
is more pronounced than that in other phases considered. 

D. Order of states at the top of the VB 

Optical and transport properties for semiconductors 
strongly depend on structure of the topmost VB. The 
order of states at the top of VB for ZnO-w is one of the 



topics which is frequently debated 0, ■ It is found that 
the normal order > Ti is obtained by LDA without SO 
coupling for all ZnX phases. Calculations within GGA 
did not changed the order. However, the order of the 
states at the top of VB for the ZnO-w is changed upon 
using the LDA+C/ approach with the semiempirical val- 
ues of the parameter U, while there appears no change 
for the other ZnX phases. For ZnO with SO coupling, 
the order of states is r7 > Tg > Tj, which is referred to as 
anomalous order, resulting from a negative SO splitting 

0- 

The variation of the structure at the top of the VB on 
U with and without SO coupling is systematically studied 
for the ZnO-z- and -w-type phases. It is found that at 
around the above considered values of U (U w 9.0 eV 
for ZnO-w and U « 8.0 eV for ZnO-z), the LDA+C/ 
interchanges the sequence of the VB states from T 7 > 

r 9 > r 7 to r 7 > r 7 > r 9 . 

Since in theoretical calculations one can freeze interac- 
tion between the valence electrons and Zn-3c? electrons by 
including the latter into the core, we studied this partic- 
ular case for ZnO-z (with SO coupling) and -w (without 
SO coupling). It is found that the the Zn-3<i electrons are 
responsible for the order of the states. On comparing the 
top of the VB structures calculated within the LDA and 
GGA approaches it is found that only small quantitative 
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FIG. 2: (a) CF splitting, SO coupling energy for ZnO-w, (b) 
SO coupling energy for ZnO-z, and (c) splitting of the states at 
the top of VB a function of U for ZnO-w. In panels a and c the 
solid and dotted lines represent calculated and experimental 
data, respectively. 
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FIG. 3: Total density of states for ZnO-z calculated from the 
LDA (solid line) and LDA+f/ (dotted line) approaches. Fermi 
energy (Ef) is set at zero of energy. 



changes have occurred. Hence, inhomogeneities in the 
electron gas do not affect the order of the states at the 
top of VB and only slightly changes the band dispersion. 

E. Effective masses 

The effective masses are calculated along T — ► A(||), 
and r — » M(-L) with and without the SO couplings (Ta- 
ble [HJ. According to the conventional notations, carrier 
masses for the ZnAT-z phases are distinguished by the in- 
dices e (electron), hh (heavy- hole), Ih (light-hole), and 
sh (split-off hole). The carrier masses for the ZnX-w 
phases are distinguished by the indices e, A, B, and C. 
The calculated CB masses m e for the ZnX-z phases are 
more isotropic than those for the ZnA"-w phases. The 
numerical values of m e for ZnO-w, ZnS-w, ZnSe-z, and 
ZnTe-z obtained by the LDA is underestimated by about 
50 % compared to experimental findings [Til u3 . Il3| . 
while those for the other ZnX phases agree fairly well 
with experimental data. GGA and LDA+J7 calculations 
only slightly improved the LDA-derived m e values for all 
ZnA" phases except ZnO, whereas for ZnO the m e values 
calculated from LDA+J7 gives better agreement with ex- 
periment. 

The heavy holes along all directions (see Table ITTfl and 
light holes along the V —> A(\\) direction are much heav- 
ier than other holes and electrons. Hence, the carrier 
transport in ZnAT is dominated by electrons, while that 
by holes can in practice be ruled out. This in turn can 
be the reason for the experimentally established large 
disparity between electron and hole mobilities I l3j . and 
also for the large optical nonlinearity in ZnO 9J. The 
hole effective masses are more anisotropic than those of 
electrons. On comparison of the values for m e in Ta- 
ble [H] one sees that the influence of SO coupling on m e 
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is very important for ZnSe-z and ZnTe-z, while for the 
other phases its effect is small. 

IV. CONCLUSION 

Electronic structure and band parameters for ZnX-w 
and -z phases are studied by first-principles calculations 
within the LDA, GGA, and LDA+J7 approaches. It is 
found that LDA underestimates the band gaps, the ac- 
tual positions of the energy levels of the Zn-3<i states as 
well as splitting energies between the states at the top 
of VB, but overestimates the crystal- field splitting en- 
ergy. The spin-orbit coupling energy is overestimated for 
ZnO-w, underestimated for ZnS-w, and comes more or 
less accurate for ZnS-z, ZnSe-z, and ZnTe-z. 

The LDA+C/ calculation as a function of U and J has 
been used to adjust the Zn-3d band derived from LDA 
to the experimentally established location of the Zn-3d 
levels determined from X-ray emission spectra. Using the 
U values corresponding to the experimentally found Ed 
location the calculated band gaps and band parameters 
are improved compared to pure LDA approach. 

The order of the states at the top of VB is systemati- 
cally examined for ZnO-z and -w phases. It is shown that 
without SO coupling the band order is normal, but it be- 
comes anomalous and Aso goes negative upon inclusion 
of SO coupling. It is found that in the LDA+U calcula- 
tions the anomalous order is maintained until U ~ 8.0 eV 
for ZnO-z and U ~ 9.0 eV for ZnO-w. Above these U val- 
ues, the band order is inverted and becomes normal and 
Acf for ZnO-w goes from positive to negative, whereas 
Aso converts to a complex quantity. It indicates that 
either the quasicubic model, within which Acf and Aso 
are calculated, do not work for this particular case or 
such order is unphysical. 

Upon excluding the interaction between the Zn-3<i and 



other valence electrons by including the former in the 
core, the order becomes anomalous. Based on these anal- 
yses one can conclude that the Zn-3<i electrons arc re- 
sponsible for the anomalous order of the states at the 
top of VB in ZnO. 

Effective masses of electrons at the conduction-band 
minimum and of holes at the valence-band maximum 
have been calculated. The heavy holes in the VB are 
found to be much heavier than the CB electrons. The 
calculations, moreover, indicate that effective masses of 
the holes are much more anisotropic than those of the 
electrons. CB electron masses for ZnO-w, ZnS-w, ZnSe- 
z, and ZnTe-z calculated within LDA arc underestimated 
by about 50 % compared with experimental data, while 
those for the other ZnX phases considered agree with 
experimental data. 

The GGA approach does not remedy the LDA derived 
error in energy gaps and band parameters. We found that 
SO coupling is important for calculation of the parame- 
ters for the z-ZnSe and z-ZnTc phases, while it is not sig- 
nificant for other ZnX phases. However, the calculated 
A CF values within the different approaches do not differ 
much except for ZnO emphasizing that Coulomb corre- 
lation effects are more pronounced for ZnO than ZnS, 
ZnSe, and ZnTe. 
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TABLE II: Effective masses of electrons and holes (in units of the free-electron mass mo) for ZnX-z and -w phases calculated 
within LDA, GGA, and LDA+C7 approaches. The results are compared to calculated and experimentally determined data cited 
in Ref. [lj| and those calculated by FP-LMTO Q, LCAO 0] and determined experimentally 0. For ZnO-w labelling of the 
effective masses was not changed with changing the order of the top VB states. 
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